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Abstract. We present a new model for the global tro- 
pospheric chemistry of inorganic bromine (Br y ) coupled 
to oxidant-aerosol chemistry in the GEOS-Chem chemical 
transport model (CTM). Sources of tropospheric Br y include 
debromination of sea-salt aerosol, photolysis and oxidation 
of short-lived bromocarbons, and transport from the strato- 
sphere. Comparison to a GOME-2 satellite climatology of 
tropospheric BrO columns shows that the model can repro- 
duce the observed increase of BrO with latitude, the northern 
mid-latitudes maximum in winter, and the Arctic maximum 
in spring. This successful simulation is contingent on the 
HOBr + HBr reaction taking place in aqueous aerosols and 
ice clouds. Bromine chemistry in the model decreases tro- 
pospheric ozone mixing ratios by < 1-8 nmol moL 1 (6.5 % 
globally), with the largest effects in the northern extratropics 
in spring. The global mean tropospheric OH concentration 
decreases by 4 %. Inclusion of bromine chemistry improves 
the ability of global models (GEOS-Chem and p-TOMCAT) 
to simulate observed 19th-century ozone and its seasonality. 
Bromine effects on tropospheric ozone are comparable in the 
present-day and pre -industrial atmospheres so that estimates 
of anthropogenic radiative forcing are minimally affected. Br 
atom concentrations are 40 % higher in the pre -industrial at- 
mosphere due to lower ozone, which would decrease by a 


factor of 2 the atmospheric lifetime of elemental mercury 
against oxidation by Br. This suggests that historical anthro- 
pogenic mercury emissions may have mostly deposited to 
northern mid-latitudes, enriching the corresponding surface 
reservoirs. The persistent rise in background surface ozone 
at northern mid-latitudes during the past decades could pos- 
sibly contribute to the observations of elevated mercury in 
subsurface waters of the North Atlantic. 


1 Introduction 

Bromine radicals (BrO x = Br+ BrO) are well-known cata- 

lysts for ozone destruction in the stratosphere (Clerbaux et 

ah, 2007), and also drive ozone depletion and mercury ox- 

idation in polar surface air (Simpson et ah, 2007; Steffen 

et ah, 2008). In the stratosphere, bromine radicals originate 

from photolysis and oxidation of halons and shorter-lived 
bromocarbons (Clerbaux et ah, 2007; Law et ah, 2007). In 
polar surface air, they originate from photochemistry involv- 

ing sea salt in aerosol or on sea ice (Simpson et ah, 2007; 

Yang et ah, 2008). There is growing evidence that bromine 
radical chemistry could be important in the background tro- 
posphere as well, with remote sensing observations typically 
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reporting 0.5-2 pmol mol - 1 BrO (Harder et al., 1998; Friel.i 
et al., 1999; Fitzenberger et al., 2000; Van Roozendael et al., 
2002; Richter et al., 2002; Platt and Honninger, 2003; Hen- 
drick et al., 2007; Theys et al., 2007, 201 1) and aircraft obser- 
vations reporting several pmol mol -1 BrO in the free tropo- 
sphere (Neuman et al., 2010; Salawitch et al., 2010). Chemi- 
cal transport model (CTM) studies show that ~ 1 pmol mol - 1 
BrO in the background troposphere would drive significant 
decrease of ozone (von Glasow et al., 2004; Lary, 2005; Yang 
et al., 2005) and could account for most of the global oxida- 
tion of elemental mercury (Holmes et al., 2006, 2010). 

Bromine radical chemistry in the troposphere is initiated 
by the production of gas-phase inorganic bromine (Br y ) from 
photolysis and oxidation of short-lived bromocarbons and 
from debromination of sea-salt aerosol (Yang et al., 2005; 
Kerkweg et al., 2008). Br y cycles between BrO x and non- 
radical reservoirs (principally HBr, HOBr, BrN 03 , BrN0 2 , 
Br 2 ) and is eventually lost by deposition. Sea salt aerosol 
debromination represents the largest source of Br y to the tro- 
posphere, with estimates between 1000 and 6000 Gg Bra -1 
constrained by the observed Br depletion in sea-salt aerosol 
relative to seawater composition (Sander et al., 2003). This 
source is mainly in the marine boundary layer where Br y 
has a short lifetime against deposition. Bromocarbons emit- 
ted by the marine biosphere including CHBn, CH 2 Br 2 , and 
CH 3 Br can release Br y in the free troposphere where its life- 
time against deposition is much longer. Global emission esti- 
mates are in the range 370-1400 Gg Bra -1 for CHBr 3 and 
57-280 Gg Bra -1 for CH 2 Br 2 (Quack and Wallace, 2003; 
Warwick et al., 2006; Law et al., 2007; Liang et al., 2010; 
Pyle et al., 2011). CH 3 Br also has a large anthropogenic 
source as an agricultural pesticide (Clerbaux et al., 2007) and 
its contribution to tropospheric Br y production has been es- 
timated at 60-80GgBra -1 (Yvon-Lewis et al., 2009). Im- 
plementation of these sources in the p-TOMCAT CTM to- 
gether with relatively well-established gas-phase Br y chem- 
istry yields a background tropospheric BrO mixing ratio of 
~ 0.5 pmol mol -1 , at the low end of observations (Yang et 
al., 2005). 

Here we incorporate a detailed simulation capability for 
tropospheric bromine coupled to oxidant-aerosol chemistry 
in the GEOS-Chem CTM (Bey et al., 2001). We evaluate this 
simulation with a recent space-based climatology of tropo- 
spheric BrO by Theys et al. (2011). We discuss the effects on 
tropospheric ozone and mercury budgets for the present-day 
and pre -industrial atmospheres. Several studies have pointed 
out the inability of current CTMs to reproduce the low sur- 
face ozone observations from the turn of the 20th century, 
suggesting that CTM-based estimates of the anthropogenic 
radiative forcing from tropospheric ozone are too low (Wang 
and Jacob, 1998; Kiehl et al., 1999; Mickley et al., 2001; 
Shindell et al., 2003; Lamarque et al., 2005; Horowitz, 2006). 
As we will show, inclusion of bromine chemistry helps to 
correct this apparent model deficiency. 
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Table 1. Global sources of tropospheric bromine in GEOS-Chem. 


Source 

Emission, 
Gg Bra -1 

Lifetime, 

d* 

Sea salt debromination 

1420 

-0 

CHBr 3 

407 

21 

CH 2 Br 2 

57 

91 

CH 3 Br 

56 

402 

Transport from stratosphere 

36 

-0 


* Global mean tropospheric lifetimes of bromocarbons determining the 
release of Bry . Values of ~ 0 are entered for the sea salt debromination 
and stratospheric sources because these sources are as Br y . 


2 GEOS-Chem bromine simulation 

The GEOS-Chem CTM (http://www.geos-chem.org) has 
been used extensively for studies of global oxidant-aerosol 
chemistry. A recent description of the model radical chem- 
istry is given in Mao et al. (2010), and a recent evaluation 
of the global tropospheric ozone simulation with sonde and 
satellite observations is given in Zhang et al. (2010). We de- 
veloped a tropospheric bromine simulation capability fully 
coupled to the standard oxidant-aerosol chemistry mecha- 
nism in version 8-02-02 of the model. The simulation in- 
cludes ten bromine species transported in the model: Br 2 , 
Br, BrO, HBr, HOBr, BrN0 2 , BrN0 3 , CHBr 3 , CH 2 Br 2 , and 
CH 3 Br. Here we describe the relevant emissions, chemistry, 
and deposition. All GEOS-Chem simulations presented here 
have horizontal resolution of 4° latitude x 5° longitude, 47 
vertical layers extending from the surface to ~ 80 km in al- 
titude, and an external chemical time step of 1 h. Transport 
is driven by GEOS-5 assimilated meteorological fields for 
2007 from the NASA Global Modeling and Assimilation Of- 
fice (GMAO). The chemical mechanism is integrated with 
the SMVGEAR solver for all tropospheric gridboxes (Ja- 
cobson and Turco, 1994; Bey et al., 2001). Linear chem- 
istry with climatological rates is used for the stratosphere 
(Bey et al., 2001), including improved treatment of ozone 
chemistry with the Linoz parameterization (McLinden et al., 
2000). Stratospheric Br y concentrations are specified from a 
model climatology as described below. 

2.1 Tropospheric sources of Br y 

Tropospheric Br y is produced in the model by debromination 
of sea salt aerosol (SSA), photolysis of CHBr 3 , and oxida- 
tion of CHBr 3 , CH?Br 2 , and CH 3 Br by OH (Table 1). SSA 
observations indicate typically a 50% loss of bromide (Br - ) 
relative to seawater composition, implying release to the at- 
mosphere as Br y (Sander et al., 2003). This release may take 
place by a set of complicated multiphase reactions producing 
Br 2 , BrCl, or HOBr (Vogt et al., 1996; Sander et al., 1999, 
2003), all of which rapidly photolyze to release BrO x radi- 
cals. Our simplified treatment of SSA debromination follows 
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Table 2a. Tropospheric bromine chemistry in GEOS-Chem: bimolecular reactions*. 


Reaction number 

Reaction 


A, cm 3 molecule 1 s 1 

-EJR, K 

R1 

CHBr 3 + OH 

— > 3Br + products 

1.35 

X 

10- 12 

-600 

R2 

Cl I 2 Br 2 + OH — > 2Br + products 

2.00 

X 

10- 12 

-840 

R3 

CH 3 Br + OH 

-* Br + products 

2.35 

X 

10- 12 

-1300 

R4 

Br + O3 -> BrO + O2 

1.60 

X 

10- 11 

-780 

R5 

Br + CH 2 0 HBr + H0 2 + CO 

1.70 

X 

10- 11 

-800 

R6 

Br + H0 2 -» 

HBr + 0 2 

4.80 

X 

10- 12 

-310 

R7 

Br + CH3CHO HBr + CH3CO3 

1.30 

X 

10- 11 

-360 

R8 

Br + (CH 3 ) 2 CO HBr + CH 3 C(0)CH 2 00 

1.66 

X 

lO-m 

-7000 

R9 

Br + C 2 H 6 HBr + C 2 H 5 00 

2.36 

X 

10“ 10 

-6411 

R10 

Br + C 3 H 8 HBr + C3H7OO 

8.77 

X 

IQ- 10 

-4330 

Rll 

Br + B1NO3 - 

— Br 2 + NO3 

4.90 

X 

10- 11 

0 

R12 

Br + NO3 — 4 

BrO + N0 2 

1.60 

X 

10- 11 

0 

R13 

HBr + OH — » 

Br + H 2 0 

5.50 

X 

10 -12 

200 

R14 

BrO + NO -> 

Br + N0 2 

8.80 

X 

10- 12 

260 

R15 

BrO + OH — » 

Br + H0 2 

1.70 

X 

10- 11 

250 

R16 

BrO + BrO - 

► 2Br + 0 2 

2.40 

X 

10- 12 

40 

R17 

BrO + BrO - 

- Br 2 + 0 2 

2.80 

X 

10-t4 

860 

R18 

BrO + H0 2 - 

> HOBr + 0 2 

4.50 

X 

10- 12 

460 

R19 

Br 2 + OH -* 

HOBr + Br 

2.10 

X 

10- 11 

240 


* The Arrhenius expression for rate constants is k = A exp (~E a /RT). Kinetic data are from Sander et al. (2011) except for 

Reactions (R7) (Atkinson et al., 2000), (R8) (King et al., 1970), (R9) and (R10) (Seakins et al., 1992), and (R1 1) (Orlando and Tyndall, 

1996). 


the Yang et al. (2005) observation-based parameterization of 
Br - depletion factors relative to seawater for particles in the 
1-10 pm diameter range. These depletion factors are applied 
to the size-dependent SSA source function in GEOS-Chem 
(Alexander et al., 2005). Uncertainties in modeled debromi- 
nation are introduced by the parameterization of bromide de- 
pletion (Sander et al., 2003) as well as the factor of 4 uncer- 
tainty in sea salt aerosol sources (Lewis and Schwartz, 2004). 
The resulting Br y is released as Br2 uniformly through the 
depth of the marine boundary layer (MBL) diagnosed from 
the GEOS-5 meteorological data. Though SSA debromina- 
tion is the largest source of tropospheric Br y in GEOS-Chem, 
removal by deposition in the MBL is fast. We find that this 
source contributes 48 % of Br y in the global free troposphere, 
the rest originating from bromocarbons. 

CHBr3 and CH?Br 2 are emitted from oceanic macroalgae 
and phytoplankton (Quack and Wallace, 2003; Yokouchi et 
al., 2005; Butler et al., 2007), and are thought to be the domi- 
nant bromocarbon precursors for tropospheric Br y (Warwick 
et al., 2006; Law et al., 2007). We use emission estimates 
from Liang et al. (2010) for different latitudinal bands and 
including enhancements in coastal regions (Carpenter et al., 
1999), based on ship cruise data (Quack and Wallace, 2003; 
Warwick et al., 2006) and refined with airborne observations 
(Liang et al., 2010). The Liang et al. (2010) emission esti- 
mates are aseasonal, but we add seasonality to the CHBr3 
source in the northern extratropics (> 30° N) as described be- 
low to better match observations. Both CHBr3 and CH 2 Br 2 


in the model are removed by oxidation by OH, and CHBr3 
is in addition removed by photolysis (Table 2). We compute 
mean tropospheric lifetimes for CHBi'3 and CH2Br2 of 21 
days and 91 days respectively (Table 1). Photochemical loss 
of CHBr3 and CH2Br2 is assumed to release all available 
Br instantaneously in our model. This approximation was 
shown to be reasonable in a modeling study by Hossaini et 
al. (2010). 

We evaluated our simulations of CHBr3 and CH2Br2 with 
observations from NASA aircraft campaigns previously used 
by Liang et al. (2010) in the evaluation of their GEOS Cli- 
mate Chemistry Model (GEOS CCM) simulation. We sam- 
pled GEOS-Chem results as averages for the coherent spatial 
domains of each campaign and for the appropriate months, 
following the standard methodology used for model evalua- 
tion of other species (Bey et al., 2001; Fischer et al., 2012). 
The observations are for years other than our 2007 simula- 
tion year but we regard interannual variability as only a small 
source of error. 

Figure 1 shows composite vertical profiles of CHBr3 and 
CH2Br2 mixing ratios for the tropics and northern extratrop- 
ics in April-June. The shapes of the vertical profiles test the 
model vertical transport modulated by atmospheric lifetime, 
including in particular the delivery to the free troposphere 
where the lifetime of Br y against deposition is long. We 
see from Fig. 1 that observed vertical gradients are steeper 
for CHBr3 than for CH2Br2, and steeper in the extratropics 
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Table 2b. Tropospheric bromine chemistry in GEOS-Chem: three-body reactions 3 . 


Reaction number 

Reaction 


Rate coefficients' 3 


R20 

Br + NO2 + M -^> 

BrN02 + M 

k 0 = 4.2 x 10“ 31 (r/300)“ 2 - 4 ; koo 

= 2.7 x 10“", F c = 0.6 

R21 

BrO + NO2 + M 

-> B1NO3 + M 

k 0 = 5.2 x 10 -31 (r/300) -3 - 2 ; k 0 0 

= 6.9 x 1(T 12 , F c = 0.6 


a The rate constant is calculated as k = {[/tq [M] / ( 1 -1- k$ iMj/^oc)] x F' 1 [ . where n = 1 1 + ( log | f j (k q | and [M] is the air number density 

(molecules cm - ''). 

k from Sander et al. (2011); and k 0 G have units of cn/ 1 molecule - " s — ' and cm" molecule - ' s - 1 respectively, and T has unit of K. 


Table 2c. Tropospheric bromine chemistry in GEOS-Chem: three- 
body reactions 3 . 


Northern Extratropics 


Tropics 


Reaction Reaction 
number 


J, s' 


-1 b 


1.1 x 1(T 6 

4.2 x ltr 2 

4.0 x 10“ 2 

2.3 x 1(T 3 

1.2 x 10“ 3 

2.2 x 1(T 4 

1.2 x 1(T 3 


R22 CHBr3 + hv — > 3Br + products 

R23 Br2 + hv — 2Br 

R24 BrO + hv -4- Br + O3 

R25 HOBr + h v -+ Br + OH 

R26 BrN02 + hv — » Br + NO2 

R27 BrN0 3 + hv —>■ BrO + N0 2 

R28 BrN03 + /iv — >■ Br + NO3 


a photolysis rates are calculated locally with the Fast-J radiative transfer model 
(Wild et al., 2000) using absorption cross-sections from Sander et al. (2011) 
except for Reaction (R26) (Scheffler et al., 1997). All quantum yields are 
assumed to be 100 % (Sander et al., 201 1). 

b global annual tropospheric mean photolysis frequencies calculated in 
GEOS-Chem. 


than in the tropics, reflecting the differences in lifetimes. The 
model reproduces well these observations. 

Figure 2 compares seasonal tropospheric columns of 
CHBr 3 and CFFB^ in the tropics and northern extratropics 
for the ensemble of NASA aircraft campaigns. This provides 
a test of emissions in the model. We find that the aseasonal 
emissions of CHBr 3 from Liang et al. (2010) overestimate 
CHBr 3 columns in the northern extratropics during winter 
and spring, and therefore we apply a seasonal correction fac- 
tor to these emissions as shown in the figure. No systematic 
bias is otherwise apparent for CHBn in the tropics or for 
CFLBr 2 . Our global CHBr 3 source is 407 Gg Bra -1 glob- 
ally (compared to 425 Gg Bra -1 in Liang et al., 2010). Pyle 
et al. (2011) recently estimated a global CHBr 3 source of 
380 Gg Bra~ constrained with observations from the OP3 
campaign over the maritime continent. Our global CFLB^ 
source is 57 Gg Br a ~ 1 , same as Liang et al. (2010). 

CHrBr has a sufficiently long lifetime to be well- 
mixed in the troposphere. Surface air observations from 
the NOAA/GMD global network show a near-uniform 
mixing ratio of 6-9 pmol mol (Montzka et al., 2003) 
and we use these as lower boundary condition in the 
model: 6.7 pmol moL 1 in the Southern Hemisphere and 
8.5 pmol mol -1 in the Northern Hemisphere. The implied 



0 0.5 1.0 1.5 2.0 2.5 0 0.5 1.0 1.5 2.0 2.5 

mixing ratio (pmol moL 1 ) 

Fig. 1. Mean vertical profiles of CHBr3 and CH2Br2 concentra- 
tions in the northern extratropics (latitudes > 30° N ) and tropics 
(30° S-30° N) for April-June. GEOS-Chem model results (in red) 
are compared to observations from NASA aircraft campaigns (in 
black) including TRACE-P, INTEX-B, and ARCTAS. See Liang et 
al. (2010) for references on the aircraft observations. Vertical pro- 
files have been averaged over coherent regions for individual mis- 
sions and then averaged again over the latitudinal domain. Horizon- 
tal bars are standard deviations in the original data. 


tropospheric source of Br y to balance the computed CH 3 Br 
sink in GEOS-Chem is 56GgBra _1 , slightly lower than 
a previous estimate of 63-78 GgBra~' (Saltzman et al., 
2004). 

Halons and other bromocarbons photolyzed or oxidized 
in the stratosphere represent an additional source of Br y to 
the troposphere. This is treated here as an upper boundary 
condition above the local model tropopause using archived 
Br y concentrations from the Liang et al. (2010) stratospheric 
simulation with the GEOS-5 CCM. The mean stratospheric 
Br y mixing ratio in that simulation is 22 pmol mol 1 . We 
use monthly mean 3-D concentration data from Liang et 
al. (2010) for individual Br y species, separately for day- 
time and nighttime. The resulting cross-tropopause source 
of tropospheric Br y in GEOS-Chem is 36GgBra _1 , as de- 
termined by difference with a sensitivity simulation using 
a zero-concentration upper boundary condition for all Br y 
species above the tropopause. 
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Table 2d. Tropospheric bromine chemistry in GEOS-Chem: heterogeneous processes. 


Reaction number 

Reaction 

Reactive uptake probability (y) a 

R29 

R30 

H5O (aerosol, cloud) 

BrN0 3 4 HOBr + HNO3 

aqueous aerosol, ice 

HOBr + HBr — » Br 2 + H 2 0 

0.3 (sea salt, liquid cloud) b , 0.8 (sulfate) c 
0.2 (sea salt, sulfate), 0.1 (ice) c 


a See Eq. (1) in text; b IUPAC evaluation (Atkinson et at., 2000b); c Sander et al. (2011); for Reaction (R30) the reactive uptake probability is 
applied to the limiting reactant. 


Table 2e. Tropospheric bromine chemistry in GEOS-Chem: 
Henry’s law constants 3 . 


Species 

Kh (298), M atm -1 

-A H/R, K 

Reference 

HBr 

0.75 (7.5 x 10 13 ) b 

10200 

Schweitzer et al. (2000) 

HOBr 

6100° 

6010 d 

Frenzel et al. (1998) 

Br b 

0.76 

4180 

Dean (1992) 

Br 

1.2 

1800 

Sander (1999) 

BrO f 

0.71 


Frenzel et al. (1998) 

BrNCh 

0.3 


Frenzel et al. (1998) 

BrNOf 

OO 


Sander (1999) 


a The Henry’s law constant is calculated as 

Kn(T) = (298) exp [— 7 ^ (t — 298 )]’ w ^ ere ^h( 298) is the value at 298 K, 

AH is the enthalpy of solution, R is the gas constant, and T is the temperature (K). 
b The Henry’s Law constant is from Schweitzer et al. (2000); The value in 
parentheses is the effective Henry’s law constant at pH 5 accounting for 

HBr( a q)/Br~ acid dissociation: K^(T) = Kfl(T) ^1 + where 

K a = 1 x 10 9 M is the acid dissociation constant (Arnaud, 1966). 
c HOBr(aq)/BrO - acid dissociation is negligible under atmospheric conditions 
(K a = 1.6 x 10 -9 M, Haag and Hoigne, 1983). 
d McGrath and Rowland (1994) 

e Bi' 2 (aq) dissociation by hydrolysis can be neglected (Beckwith et al., 1996). 
f Estimated by analogy with CIO (Sander et al., 201 1). 
g B 1 NO 3 uptake by liquid clouds is modeled by a reactive uptake probability 
(Table 2d) and is not limited by Henry’s law solubility. For modeling dry deposition 
an infinite Henry’s law constant is assumed. 


Other tropospheric Br y sources are neglected as unim- 
portant on a global scale. Kerkweg et al. (2008) estimated 
that photo-oxidation of bromocarbons other than CHBrq, 
CH 2 Br 2 , and CH(Br accounts for 24 Gg Bra' 1 of tropo- 
spheric Br y . Pyle and Mather (2009) estimated a global vol- 
canic source of 5-15 Gg Bra -1 . Photochemical bromine re- 
lease from sea ice and salt beds can lead to high local concen- 
trations of BrO in surface air (Stutz et al., 2002; Simpson et 
al., 2007) but the volumes affected are negligible on a global 
scale. 

2.2 Chemical cycling of Br y 

Table 2 lists the chemical mechanism for tropospheric 
bromine in GEOS-Chem. Rate constants, heterogeneous re- 
action coefficients, and photolysis cross-sections are taken 
from Sander et al. (2011) unless otherwise specified. Pho- 
tolysis rates are calculated online with the Fast-J radiative 
transfer model (Wild et al., 2000). Heterogeneous bromine 
chemistry in aqueous aerosols includes hydrolysis of BrNOa 



J FMAMJ JASONDJ FMAMJ J ASOND 


Month 

Fig. 2. Comparison of simulated (red) and observed (black) tropo- 
spheric columns (molecules cm -2 ) of CmBm and CHBr3, for the 
northern extratropics (> 30° N) and the tropics (30° S-30° N). Col- 
umn data are integrals of mean vertical profiles up to 10 km alti- 
tude from a compilation of NASA campaigns including STRAT, 
TRACE-P, INTEX, AVE, TC4, and ARCTAS in the northern extra- 
tropics and STRAT, PEM-Tropics, TRACE-P, INTEX, Pre-AVE, 
AVE, and TC4 in the tropics. References for the observations are 
in Liang et al. (2010). Vertical bars are the standard deviations 
in the original data. The dashed red line shows the model results 
for CHBr3 in the northern extratropics using the original aseasonal 
emissions from Liang et al. (2010); in order to match the observa- 
tions in GEOS-Chem we applied seasonal scaling factors to CHB13 
emissions in that region (solid red line). 

(Reaction R29) and reaction between HOBr and HBr (Reac- 
tion R30). Both are treated with a standard reactive uptake 
probability (y) formulation, resulting in an effective gas- 
phase loss rate constant to a monodisperse aerosol (radius 
a) given by 



Here D is the molecular diffusion coefficient in air, v 
is the mean molecular speed of the gas, and A is the 
aerosol surface area concentration per unit volume of air 
(cm 2 cm -3 ) (Schwartz, 1986; Jacob, 2000). GEOS-Chem 
simulates explicitly the mass concentrations of different 
aerosol types, and we integrate k over the prescribed rela- 
tive humidity-dependent size distributions as described by 
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Martin et al. (2002). Reactive uptake probabilities for both 
Reactions (R29) and (R30) are based on existing laboratory 
evidence (Table 2d) but with limited information on the sen- 
sitivity to temperature and aerosol composition. For BrNC>3 
hydrolysis Reaction (R29) on liquid cloud droplets, the cloud 
surface area concentration A is calculated from GEOS-5 liq- 
uid water content data, assuming effective droplet radii of 
10 pm and 6 pm for marine and continental clouds respec- 
tively (Park et al., 2004; Fu et al., 2008). For HBr+HOBr 
(Reaction R30) on ice cloud surfaces, we use A — 2 x 
10 ~ 4 /0.9 cm 2 cm -3 (Lawrence and Crutzen, 1998) applied to 
the local ice water content I (cm 3 cm -3 ) from the GEOS-5 
data. 

Reaction (R30) is an important pathway for recycling 
bromine radicals in the troposphere (von Glasow et al., 
2004; Yang et al., 2010). Although the kinetics are uncertain, 
there is substantial evidence for this reaction. Abbatt (1995) 
measured y > 0.25 for HBr + HOBr in sulfuric acid solu- 
tion (69wt%) with HBr present in excess. Laboratory mea- 
surements show that HOBr reactive uptake on deliquescent 
NaBr and NaCl aerosol has y > 0.2 for pH < 7 (Abbatt and 
Waschewsky, 1998; Fickert et al., 1999; Wachsmuth et al., 
2002). Fickert et al. (1999) determined the relative produc- 
tion of Br2 vs. BrCl for HOBr reacting with aqueous salt 
solutions containing both Br - and Cl - , and found >90% 
Br 2(g) production for [Br - ]/[Cl - ] ratios typical of seawa- 
ter. Mochida et al. (1998) measured y — 0.18 ±0.02 for 
HOBr reactive uptake on solid KBr. On ice surfaces with 
excess HBr, Abbatt (1994) measured y— 0.12 ±0.03 for 
HOBr + HBr at 228 K, and Chaix et al. (2000) measured 
y = 0.15 ±0.03 at 205 K. Based on this ensemble of evi- 
dence and the recommendation of Sander et al. (2011), we 
assume that Reaction (R30) proceeds with y = 0.2 for sul- 
fate and sea salt aerosol (presumed aqueous) and y = 0.1 for 
ice surfaces, where y is applied to the locally limiting reac- 
tant. This simplified approach assumes that the non-limiting 
gas is at equilibrium with the aerosol phase. 

2.3 Sinks of Br y 

Individual Br y species are removed by wet and dry deposi- 
tion. Wet deposition in GEOS-Chem takes place by scaveng- 
ing in convective updrafts and by in- and below-cloud scav- 
enging from large-scale and convective precipitation (Liu et 
al., 2001). A standard test for wet deposition in global CTMs 
is simulation of 210 Pb aerosol, for which the source and 
atmospheric concentrations are relatively well constrained 
(Balkanski et al., 1993). Our version of GEOS-Chem yields a 
tropospheric lifetime of 210 Pb against deposition of 9.5 days 
(L. T. Murray, personal communication, 2011), consistent 
with a best observational estimate of 9 days (Liu et al., 2001). 

Different gas scavenging efficiencies are used in GEOS- 
Chem for warm liquid clouds ( T > 268K), mixed clouds 
(248 < T < 268 K), and cold ice clouds ( T < 248 K). For 
warm clouds, partitioning into the aqueous phase is deter- 
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mined by the Henry’s law constant (Table 2e) and the liq- 
uid volume fraction of the precipitating water. Precipitation 
in mixed clouds is assumed to take place by riming, and 
the scavenging efficiency depends on the retention efficiency 
when the cloud droplets freeze (Stuart and Jacobson, 2003). 
For HBr we assume a 100 % riming retention efficiency, fol- 
lowing recommendations from Stuart and Jacobson (2003) 
for gases with large effective Henry’s Law constants. For 
other Br y species we assume no retention during riming. We 
find little sensitivity to the assumed retention efficiency, con- 
sistent with previous GEOS-Chem studies for scavenging of 
ammonia (Wang et al., 2008; Fisher et al., 2010). We assume 
no scavenging from cold clouds. HBr and HOBr undergo het- 
erogeneous chemistry at the surface of ice particles as dis- 
cussed in Sect. 2.2. 

Dry deposition velocities are calculated locally in GEOS- 
Chem following a standard resistance-in-series model with 
surface resistance determined by the Henry’s law constant, 
surface type, and meteorological conditions (Wesely, 1989; 
Wang et al., 1998). The resulting global annual mean depo- 
sition velocities are 0.94 cm s -1 for B1NO3, 0.93 cm s -1 for 
HBr, and 0.40 cm s -1 for HOBr. 


3 Global budget of tropospheric Br y 

Figure 3 shows the global annual mean budget and cy- 
cling of tropospheric Br y in our simulation. The BrO mx- 
ing ratio is 0.32 pmol mol -1 , which implies a daytime mean 
of 0.64 pmol mol - 1 since nighttime concentrations are near 
zero. This is at the low end of the oft-cited (but poorly 
constrained) 0.5-2 pmol mol -1 range in the observations; 
more detailed comparison with satellite observations will 
be presented in Sect. 4, and comparison against BrO ob- 
servations taken at Cape Verde will be presented by Stone 
and Evans (2012). The global mean Br y mixing ratio is 
3.2 pmol mol -1 , with HBr and HOBr as the principal reser- 
voirs accounting respectively for 34 % and 28 % of total Br y . 
Br y has a lifetime of 7 days against deposition, with HBr 
accounting for 55 % of that sink and HOBr for 40 %. Chemi- 
cal lifetimes are short relative to deposition; that of HBr (the 
longest-lived reservoir) is only 6h. Thus the BrO x radical 
concentrations are effectively maintained by chemical recy- 
cling from non-radical reservoirs. 

We see from Fig. 3 that the HBr+HOBr heterogeneous 
Reaction (R30) plays an important role in the recycling of 
BrO x from HBr. A sensitivity simulation without this reac- 
tion indicates a factor of 2 decrease in the BrO mixing ratio 
to 0.1 5 pmol mol -1 , with HBr accounting for 70% of Br y ; 
the lifetime of Br y is shorter in that simulation (5 days) be- 
cause of the high water solubility of HBr (Table 2e). The im- 
portance of Reaction (R30) reflects the comparable concen- 
trations of HBr and HOBr (Fig. 3), as rapid HOBr production 
by the BrO + H0 2 Reaction (R18) compensates for the much 
longer lifetime of HBr (6 h) than HOBr (20 min). 
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Deposition Sea salt aerosol CHBr 3 CH 2 Br 2 CH 3 Br 
debromination 

Fig. 3. Global annual mean budget of tropospheric inorganic 
bromine (Bry) in GEOS-Chem. The main reactions are indi- 
cated. Inventories are given as masses (Gg Br), with mixing ratios 
(pmolmol - 0 in brackets. Rates are given in units of GgBrs - . 
Read 6.3(— 5) as 6.3 x 10 -5 . HBr accounts for 55 % of Br y loss by 
deposition. Sea-salt aerosol debromination is the dominant global 
source of Br y but is mainly confined to the marine boundary layer 
where Br y has a short lifetime against deposition. It accounts for 
48 % of the Br y source in the global free troposphere. 


Figure 4 shows the zonal annual mean BrO, Br, and Br y 
mixing ratios in the model. The highest concentrations of 
Br y and BrO are in the marine boundary layer at southern 
mid-latitudes due to the large source from sea-salt aerosol 
debromination. In the rest of the troposphere Br y is relatively 
uniform (2^-pmolmol -1 ), with a minimum in the tropical 
upper troposphere due to scavenging in convective updrafts. 
Mixing ratios increase with altitude in the stratosphere where 
the boundary conditions from Liang et al. (2010) are applied. 

BrO and Br concentrations in Fig. 4 show a general in- 
crease with altitude, and BrO shows in addition an increase 
with latitude. The budget of BrO is largely that of BrO x since 
[BrO] J>> [Br], From Fig. 3 we see that the dominant sink 
of BrO x is the BrO + lKL Reaction (R18), which is fastest 
in the tropical lower troposphere where HO 2 concentrations 
are highest. The dominant source of Br is photolysis of BrO 
(Fig. 3), which cancels the latitudinal dependence of BrO 
concentrations and drives the strong increase of Br concen- 
trations with altitude. 



Fig. 4. Zonal annual mean concentrations (pmol mol 1 ) of BrO, Br, 
and Br y in GEOS-Chem. 


4 Comparison to satellite observations of tropospheric 
BrO 

Figure 5 compares simulated BrO tropospheric columns to a 
GOME-2 satellite data set retrieved by Theys et al. (2011) 
with monthly resolution for 2007 and for six latitudinal 
bands. Theys et al. (2011) fitted total slant columns of BrO 
to the solar backscattered radiation spectra. They then re- 
moved the stratospheric contribution by using correlations of 
stratospheric BrO with ozone and NO 2 from the BASCOE 
3D-CTM (Theys et al., 2009) and applying these correla- 
tions to the local GOME-2 observations of ozone and strato- 
spheric NO 2 columns. They converted the tropospheric slant 
columns to vertical columns with an air mass factor assum- 
ing a Gaussian-shaped vertical distribution of BrO peaking 
at 6 km altitude. The grey shading in Fig. 5 indicates the 
total observational error estimated by Theys et al. (2011), 
including contributions from the air mass factor, surface 
albedo, cloud fractions, stratospheric BrO column removal, 
and cloud top heights. 

The GOME-2 data are for ~ 09:30 local time (LT) and 
we sample GEOS-Chem at 09:00-1 1:00 LT for comparison. 
GEOS-Chem data in Fig. 5 include both the actual tropo- 
spheric vertical columns (red) and the columns corrected for 
the Theys et al. (201 1) air mass factors using their tabulated 
scattering weights including for partly cloudy scenes. The 
correction is uncertain because the tabulation of scattering 
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weights is sparse, but the effect is relatively small except 
at southern mid-latitudes where the model BrO has a large 
contribution from the boundary layer. Also shown in Fig. 5 
are actual tropospheric vertical columns sampled at 09:00- 
1 1 :00 LT from the p-TOMCAT CTM bromine simulation of 
Yang et al. (2010), (blue) and from GEOS-Chem without 
the HBr + HOBr Reaction (R30) (dashed). p-TOMCAT and 
GEOS-Chem have completely different heritages. Specific 
differences in the p-TOMCAT simulation of bromine include 
the use of a smaller rate constant for HBr + HOBr (y =0.02 
instead of 0.2, and no reaction on cloud ice particles), and 
consideration of a blowing snow source in polar regions dur- 
ing spring. 

The standard GEOS-Chem simulation underestimates the 
GOME-2 observations by 30 % on a global mean basis, ex- 
cluding the poles (40 % including the air mass factor correc- 
tion). This is marginally significant considering the observa- 
tion error but might indicate a low bias in the simulated emis- 
sions or heterogeneous recycling of Br y . The observations 
show a positive BrO gradient with increasing latitude that 
is well reproduced in GEOS-Chem contingent on the (Re- 
action R30) HBr + HOBr reaction. Although Reaction (R30) 
is a photochemical pathway inasmuch as HOBr and HBr are 
mainly produced photochemically, it is most effective when 
radiation is weak since aerosol uptake of HOBr competes 
with HOBr photolysis. Thus the effect of Reaction (R30) 
for BrO x recycling is relatively small in the tropics but very 
large in the extratropics. It has a strong effect on seasonality 
in the extratropics, enabling GEOS-Chem to reproduce the 
observed winter maximum at northern mid-latitudes, in con- 
trast to p-TOMCAT where Reaction (R30) is less important 
due to a seasonally constant non-sea salt aerosol loading and 
lower y. 

We conducted several sensitivity simulations to examine if 
uncertainties in gas-phase bromine kinetics as estimated by 
Sander et al. (201 1) could account for the 30 % global model 
underestimate of BrO relative to the Theys et al. (2011) data. 
We found the effects to be relatively small. The largest re- 
sponse was from the Br + CH20 Reaction (R5). Decreasing 
the pre-exponential factor for that reaction by 30 % and in- 
creasing the activation energy by 30% (error standard de- 
viations from Sander et al., 2011) resulted in a 19% global 
increase in simulated BrO. 

A major feature in the BrO observations of Theys et 
al. (2011) is the strong Arctic spring maximum. This 
is reproduced by both GEOS-Chem and p-TOMCAT but 
for different reasons. In GEOS-Chem it is due mainly to 
fast HBr+HOBr recycling (Reaction R30) (yHOBr = 0.2 in 
GEOS-Chem compared to ynoBr = 0.02 in p-TOMCAT for 
aqueous aerosol), with no springtime emissions of Br y from 
sea ice. The spring shift in GEOS-Chem (relative to the win- 
ter maximum at northern mid-latitudes) is due to the need 
for some insolation to drive HOBr production. p-TOMCAT 
simulates a spring maximum with slower HBr + HOBr re- 
cycling by including a blowing-snow source for Br y (Yang 
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Fig. 5. Seasonal variation of mean tropospheric BrO columns in 
different latitudinal bands. GOME-2 observations from Theys et 
al. (2011) are compared to model values from GEOS-Chem and 
from p-TOMCAT (Yang et al., 2010). GEOS-Chem values are 
shown for the actual vertical columns (red) and for vertical columns 
corrected by the GOME-2 air mass factors (green). Also shown are 
results from a sensitivity simulation without the HBr+HOBr het- 
erogeneous Reaction (R30). p-TOMCAT values are actual vertical 
columns. The models were sampled at 09:00-1 1 :00 local time (LT), 
corresponding to the GOME-2 viewing time of ~ 09:30 LT. Grey 
shading shows the observational error on the means reported by 
Theys et al. (2011) and vertical bars indicate one standard devia- 
tion of their spatial averaging. The BrO mixing ratios on the right 
axis scale are derived from the BrO columns by dividing by the an- 
nual mean tropospheric air density columns for the corresponding 
latitude bin in GEOS-Chem. 


et al., 2008, 2010). The Arctic column amounts in GEOS- 
Chem are dominated by the free troposphere, while those in 
p-TOMCAT have a large boundary layer contribution. In-situ 
observations of BrO during the ARCTAS aircraft campaign 
over Alaska in April 2008 found significant mixing ratios in 
the free troposphere, with a mean of ~ 2 pmol mol 1 (Salaw- 
itch et al., 2010). In comparison, GEOS-Chem simulates a 
09:00-17:00 LT mean of 1.3 pmol mol -1 BrO for the Arc- 
tic free troposphere in April. Both models likely oversim- 
plify photochemistry in polar spring, as local emissions of 


Atmos. Chem. Phys., 12, 6723-6740, 2012 


www.atmos-chem-phys.net/12/6723/2012/ 


J. P. Parrella et al.: Tropospheric bromine chemistry 

Br y and the extent of heterogeneous recycling remain poorly 
constrained (Simpson et al., 2007). 

Figure 5 shows little seasonal variation of BrO column 
in the tropics, either in the observations or the models. 
The models are too low. Observations at southern mid- 
latitudes show concentrations comparable to northern mid- 
latitudes (confirming the dominant natural origin of Br y ) 
but weaker seasonal dependence as might be expected from 
lower aerosol concentrations to drive (Reaction R30). 

5 Impact on tropospheric ozone and OH 

Tropospheric bromine chemistry drives several catalytic cy- 
cles for ozone loss in the troposphere. The three most impor- 
tant are: 

Mechanism 1 : BrO + HO 2 

Br + O 3 — » BrO + O 2 (R5) 

BrO + H0 2 — * HOBr + 0 2 (R14) 

HOBr + hv — ^ Br + OH (R25) 

net: O 3 + HO 2 — > OH + 2 O 2 

with closure for the hydrogen oxide radicals (HO x = OH + 
H + HO 2 + organic peroxy and oxy radicals) provided for ex- 
ample by oxidation of CO, 

co + oh-^co 2 + ho 2 . 

Mechanism 2: BrN 03 + H 2 O 

Br + O 3 — » BrO + O 2 (R5) 

BrO + NO 2 + M — * BrN0 3 + M (R20) 

BrN0 3 — a - r ° - °i HOBr + HNO 3 (R29) 

HOBr + hv — ^ Br + OH (R25) 

net: O 3 + NO 2 + H 20 ( a q ) — > O 2 + HNO 3 + OH . 

Mechanism 3: HOBr + HBr 


Br + O 3 — > 

BrO + O 2 

(Rl) 

BrO + HO 2 - 

-> HOBr + 0 2 

(R 2 ) 

O 

Br + RH— * 

HBr + R0 2 

(R5-R10) 


aerosol or ice _ 


HOBr + HBr 

> Br 2 + H 2 0 

(R30) 

Br 2 + hv — > 

► 2Br 

(R23) 


net : O 3 + HO 2 + RH — > O 2 + H 2 O + RO 2 . 

A common step in all three mechanisms is the formation of 
HOBr. Loss of HOBr is either by photolysis or by reaction 
with HBr (Fig. 3), and both of these sinks drive catalytic 
ozone loss (Mechanisms 1 and 3). Thus Br-catalyzed ozone 
loss is limited by the rate of HOBr production. This is in 
contrast to the polar springtime boundary layer, where the 
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BrO + BrO self-reaction is the dominant catalytic ozone sink 
due to very high BrO concentrations. We see from Fig. 3 that 
the self-reaction of BrO is 50 times slower than the reac- 
tion with HO 2 on the global scale. The Br0 + H 02 Reac- 
tion (R14) accounts for 95 % of global tropospheric HOBr 
production while BrN 03 hydrolysis accounts for the rest. 
Since 92 % of HOBr is photolyzed. Mechanism 1 is respon- 
sible for about 90 % of Br-catalyzed ozone loss. However, 
Mechanism 2 has the important secondary effect of provid- 
ing a sink for NO x and thus slowing down ozone production. 

Table 3 shows the global impact of bromine chemistry 
on the tropospheric ozone budget. The budget is calculated 
following standard procedure as that of the odd oxy- 
gen family O x = O 3 + NO 2 + 2 NO 3 + peroxyacylnitrates 
+ HNO 4 + 3N 2 0 5 + HNO 3 + BrO + BrN0 2 + 2BrN0 3 
(Crutzen and Schmailzl, 1983; Wu et al., 2007) to ac- 
count for fast cycling of minor species with ozone. The 
tropospheric production rate of ozone in GEOS-Chem 
decreases by 4.1% when bromine chemistry is included, 
and the ozone lifetime decreases by 2.9%. We find a 
6.5 % decrease in the global tropospheric ozone burden. A 
previous model study by von Glasow et al. (2004) found 
a larger (10-15%) effect of bromine chemistry on the 
tropospheric ozone budget because they imposed higher 
bromine concentrations. Yang et al. (2010) found bromine to 
decrease the ozone burden by 25.5 Tg ( 8 . 6 %), close to our 
result of 24 Tg (6.5%). 

We see from Table 3 that bromine chemistry also causes 
a 4 % decrease in the global concentration of OH, the main 
tropospheric oxidant. This appears to mainly reflect the de- 
creases in ozone and NO x , and is partly compensated by 
the source of OH from Mechanism 2. Applying our global 
mean fields of OH and temperature with kinetics from Sander 
et al. ( 2011 ), we calculate a methylchloroform tropospheric 
lifetime against oxidation by OH in GEOS-Chem of 4.79 yr 
with bromine and 4.65 yr without. Though both are shorter 
than estimates from observations, 5.5 to 6 yr (Spivakovsky et 
al., 2000; Prinn et al., 2005), bromine decreases the model 
bias. Fast conversion of HO 2 to H 2 O in aerosols, not in- 
cluded in our GEOS-Chem version, could fully correct the 
bias (Mao et al., 2010, 2012). 

Figure 6 shows the seasonal zonal mean decreases in tro- 
pospheric ozone mixing ratios due to bromine chemistry. 
Concentrations decrease by < 1-8 nmol mol -1 depending on 
region and season. The impact is largest in the northern ex- 
tratropics in spring, reflecting a combination of elevated BrO 
(Fig. 5) and ozone. The general increase with latitude in 
the effect of bromine on ozone is consistent with the pre- 
vious model studies of von Glasow et al. (2004) and Yang et 
al. (2010). We find little zonal variability in the effect. 

A number of studies have previously evaluated the global 
ozone simulation in GEOS-Chem with observations from 
surface sites, sondes, and satellites (Sauvage et al., 2007; 
Nassar et al., 2009; Zhang et al., 2010). These have 
shown patterns of over- and underestimates that can vary 
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Table 3. Global tropospheric ozone (O x ) budget in GEOS-Chem a . 


Chemical production, Chemical loss, Deposition STE b Tropospheric Lifetime [OH] 

P(O x ) L(O x ) burden 

Tga -1 Tga -1 Tga -1 Tga -1 Tg days 10 6 molecules cm -3 

No bromine 5110 4480 948 349 370 27.8 1.23 

With bromine 4910 4350 879 355 346 27.0 1.18 


a Global annual mean budget for the odd oxygen (O x ) chemical family defined as O x = O3 + NO2 + 2NO3 + peroxyacylnitrates + HNO4 + 3N2O5 + HNO3 + BrO + B1NO2 + 
2BrNC>3 (Crutzen and Schmailzl, 1983 ; Wu et al., 2007 ), with molecular weight taken to be that of ozone, 
b Stratosphere-to-troposphere exchange is inferred from mass-balance: STE = L(O x ) + Deposition - P(O x ). 
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Fig. 6. Decrease in tropospheric ozone (nmol mol 1 ) from bromine 
chemistry. The figure shows zonal seasonal mean differences 
between GEOS-Chem simulations including vs. not including 
bromine chemistry. Shading indicates the stratosphere. 

significantly in successive model versions due to changes in 
model emissions, chemistry, and transport (Wu et al., 2007). 
Zhang et al. (2010) presented an extensive global comparison 
to ozonesonde and satellite measurements using a version of 
GEOS-Chem similar to ours (v8-02-03). They found that the 
model was too high in southern mid-latitudes and the north- 
ern subtropics but too low in the tropics, with negligible bias 
in the northern extratropics. Adding bromine chemistry in 
that simulation would produce mixed results in comparisons 
to observations. However, the latest version of GEOS-Chem 
(v9; http://www.geos-chem.org) includes improved isoprene 
chemistry that corrects the ozone bias in the tropics but 
drives a systematic model overestimate of 5-10 nmol mol -1 
in the northern extratropics (Paulot et al., 2009a, b). Includ- 
ing bromine chemistry in that latest version could largely cor- 
rect the discrepancy in the extratropics and this will need to 
be investigated in future work. 

6 Impact on pre-industrial ozone and radiative forcing 

A long-standing problem in global modeling of tropospheric 
ozone has been the inability of models to reproduce the very 
low ozone concentrations measured in surface air worldwide 


Montsouris, France (1876 - 1886) 



Month 

Fig. 7. Mean seasonal variation of 19th-century ozone in surface 
air at Montsouris, France (50° N, 2° E). 1876-1886 observations re- 
analyzed by Volz and Kley (1988) (black) are compared to preindus- 
trial model simulations using GEOS-Chem (red) and p-TOMCAT 
(blue), with bromine chemistry (solid) and without (dotted). 


at the turn of the 20th century (Wang and Jacob, 1998; Mick- 
ley et al., 2001; Shindell et al., 2003; Lamarque et al., 2005; 
Horowitz, 2006). Calibration of these older measurements 
is controversial (Marenco et al., 1994) but the Montsouris 
(Paris) measurements are considered reliable (Volz and Kley, 
1988). To examine the potential for bromine chemistry to ad- 
dress this discrepancy we conducted a GEOS-Chem sensitiv- 
ity simulation for the pre-industrial atmosphere removing all 
anthropogenic sources (including NO x from fertilizer use), 
and reducing methane from 1700 to 700 nmol mol -1 . We de- 
creased biomass burning to 10% of its present value (Wang 
and Jacob, 1998), though this has little impact on our re- 
sults. Bromine sources were held to present-day estimates, 
except for CH 3 Br and the stratospheric boundary condition. 
We used 5 pmol moL 1 CHBr 3 in surface air worldwide based 
on ice core records (Saltzman et al., 2004). Mixing ratios of 
stratospheric Br y species were scaled from a total of 22 to 
1 2 pmol moL Br y (Liang et al., 2010; Montzka et al., 201 1). 
Sea salt debromination rates may have increased between the 
preindustrial and present-day due to enhanced particle acid- 
ity from anthropogenic emissions (Sander et al., 2003) but 
this effect is not considered here. 
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Figure 7 compares the 1876-1886 Montsouris observa- 
tions to simulated pre-industrial ozone from GEOS-Chem 
and p-TOMCAT (Yang et al., 2010), with and without 
bromine chemistry. p-TOMCAT applies similar conditions to 
GEOS-Chem for simulating the pre-industrial atmosphere. 
The simulations without bromine chemistry overestimate 
the mean observed ozone by 6.2 and 3.1nmolmol -1 in 
GEOS-Chem and p-TOMCAT respectively. Both simula- 
tions show a seasonal maximum in late winter and early 
spring, which is a standard extratropical feature in mod- 
els of natural ozone (Wang and Jacob, 1998; Monks, 2000; 
Mickley et al., 2001) but is not seen in the Montsouris ob- 
servations. Including bromine chemistry decreases annual 
mean ozone at Montsouris by 2.6 and 1.9 nmol mol -1 in 
GEOS-Chem and p-TOMCAT, respectively, effectively cor- 
recting the overestimate in p-TOMCAT and reducing it sub- 
stantially in GEOS-Chem. The decrease is least in sum- 
mer (0.9-1 .4 nmol mol -1 ) and greatest in winter-spring (2.3- 
5.1 nmol mol -1 ), which suppresses the winter-spring max- 
imum in both models. Although there are significant dif- 
ferences in simulated ozone between GEOS-Chem and p- 
TOMCAT, the impact of bromine chemistry on both models 
is consistent. 

Saiz-Lopez et al. (2012) recently suggested that halogen 
chemistry could reduce the radiative forcing from anthro- 
pogenic tropospheric ozone since pre-industrial times by 
—0.1 W m -2 or about 30%. However, they only considered 
the impact of this chemistry on the present-day atmosphere. 
Considering that tropospheric halogens are mainly natural, 
they would also decrease ozone in the pre-industrial atmo- 
sphere. We find in GEOS-Chem that although the sources of 
bromine are similar in the present and pre-industrial atmo- 
spheres, the chemistry differs in two principal ways. First, 
lower methane in the pre-industrial atmosphere leads to 
lower CH 2 O and hence suppressed production of HBr by 
Reaction (R5). Second, lower aerosol loadings suppress the 
HBr + HOBr heterogeneous Reaction (R30). These two ef- 
fects have competing influences on tropospheric BrO and we 
find in GEOS-Chem that they largely cancel. Thus the effects 
of bromine chemistry on ozone are similar in the present-day 
and pre-industrial atmospheres. We find that the increases in 
the total tropospheric ozone burden from pre-industrial to 
present are 119 and 113Tg in the simulations without and 
with bromine chemistry, respectively. Considering that the 
radiative forcing from tropospheric ozone depends roughly 
linearly on the change in burden (Mickley et al., 1999), we 
conclude that bromine chemistry decreases ozone radiative 
forcing by about 5 %. 

The impact of bromine chemistry on pre-industrial US sur- 
face ozone is of interest because this natural background 
represents a policy reference point for assessing the risk 
associated with anthropogenic emissions (US Environmen- 
tal Protection Agency, 2006). Zhang et al. (2010) recently 
reported a mean natural ozone background in US surface 
air of 18 ±6 nmol mol -1 using a high-resolution version of 


GEOS-Chem. We find in our simulation that bromine chem- 
istry decreases natural surface ozone over the US from 16 
to 13 nmol mol -1 on an annual mean basis, i.e., an effect of 
3 nmol mol -1 or 19 %. 

7 Implications for the atmospheric lifetime of mercury 

Mercury in the atmosphere is mostly emitted as gaseous 
Hg(0), and is removed by oxidation followed by rapid de- 
position of the water-soluble Hg(II) compounds (Lindberg et 
al., 2007). The observed atmospheric variability of Hg(0) im- 
plies a lifetime against deposition of 0.5-1 yr, with a possibly 
shorter lifetime against oxidation depending on whether at- 
mospheric Hg(II) can be reduced back to Hg(0) (Slemr et al., 
1981, 1985; Bergan et al., 1999; Holmes et al., 2010). Most 
models of atmospheric mercury have assumed that oxidation 
of Hg(0) is driven by OH and ozone, but these reactions are 
now thought to be too slow to be of atmospheric relevance 
(Calvert and Lindberg, 2005; Hynes et al., 2009). Goodsite 
et al. (2004) showed that rapid oxidation of Hg(0) to Hg(II) 
by Br atoms could take place by 

Hg + Br + M — s- HgBr + M (R31) 

HgBr 4- Hg + Br (R32) 

HgBr + X 4- HgBrX , (R33) 

where X is a radical, such as OH or Br, that converts unsta- 
ble Hg(I) to stable Hg(II) in a 3-body reaction. Under tro- 
pospheric conditions. Reaction (R31) operates in the low- 
pressure limit whereas Reaction (R33) operates in the high- 
pressure limit (Goodsite et al., 2004, 2012; Donohoue et al., 
2006). Using the kinetic data from Goodsite et al. (2004) 
and Donohue et al. (2006), together with Br concentration 
fields from the p-TOMCAT model (Yang et al., 2005) and 
OH fields from GEOS-Chem, Holmes et al. (2006, 2010) 
showed that oxidation by Br atoms could represent the dom- 
inant global atmospheric sink for Hg(0). Here we revisit this 
calculation using our computed Br concentrations and draw 
implications for possible changes in the Hg(0) atmospheric 
lifetime between pre-industrial and present conditions. 

The lifetime r of Hg(0) against oxidation to Hg(II) by Re- 
actions (R3 1)-(R33) can be expressed assuming HgBr steady 
state by Holmes et al. (2006): 

T _ £32 + fc 33 ,Br[Br] + k33 : on[OH] 

^31 [Br] [M] (k 33 ,Br[Br] + k 3 3 ,OH[OH]) ’ 

where we have assumed that X in Reaction (R33) can be ei- 
ther Br or OH, following Holmes et al. (2006, 2010). Using 
the kinetic data of Holmes et al. (2010) with updates from 
Goodsite et al. (2012), and assuming a uniform tropospheric 
Hg(0) concentration, we derive from our GEOS-Chem model 
fields a global tropospheric Hg(0) lifetime of 0.58 yr against 
oxidation by Br. In comparison. Holmes et al. (2010) calcu- 
lated a lifetime of 0.50 yr in a GEOS-Chem Hg(0) simulation 
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using the Yang et al. (2005) Br mixing ratios in the free tropo- 
sphere and assuming 1 pmol mol -1 BrO in the marine bound- 
ary layer. We repeated our calculation using their archived 
Br and OH fields and obtain a lifetime of 0.41 yr with the 
Holmes et al. (2010) kinetics, and 0.65 yr with the kinetics 
update from Goodsite et al. (2012). Our estimate of 0.41 yr 
can be viewed as consistent with the Holmes et al. (2010) 
estimate of 0.50 yr in view of the uncertainty involved with 
assuming a uniform Hg(0) concentration. Our 1 1 % shorter 
lifetime (0.58 vs. 0.65 yr) reflects a 13 % greater mean tropo- 
spheric Br concentration in our GEOS-Chem simulation. 

The mean tropospheric Br atom concentration in our pre- 
industrial simulation is 38 % greater than present-day due to 
lower ozone suppressing the Br + O3 Reaction (R4). This in- 
crease in Br causes the Hg(0) lifetime against oxidation to 
decrease to 0.35 yr, 40 % lower than present-day. Our work 
thus suggests that the increase in ozone since pre-industrial 
times may have increased the atmospheric residence time of 
Hg(0) by 66 %. With an atmospheric lifetime of four months, 
mercury in pre-industrial times would be less of a global pol- 
lutant than today. Most of the observed increase in ozone has 
taken place since 1950 (Marenco et al., 1994), by which time 
there had already been large anthropogenic mercury emis- 
sions at northern mid-latitudes from mining, industry, and 
combustion (Streets et al., 2011). 

Our work highlights a potentially important link between 
ozone trends and the lifetime of Hg(0) through partitioning 
of the Br y family. This result is contingent on Br playing a 
dominant role in Hg(0) oxidation. Changes in sea salt de- 
bromination between the preindustrial and present as well as 
our simplified treatment of heterogeneous bromine chemistry 
represent additional uncertainties that might affect our Hg(0) 
lifetime estimates. 


8 Conclusions 

We have developed a simulation of tropospheric bromine 
coupled to oxidant-aerosol chemistry in the GEOS-Chem 
global chemical transport model (CTM). Our goal was to 
better understand the effect of bromine chemistry on tro- 
pospheric ozone and mercury, for both the present and pre- 
industrial atmospheres. 

We found that standard gas-phase mechanisms for 
bromine chemistry are unable to reproduce observed levels 
of tropospheric BrO, including the latitudinal and seasonal 
variations reported in the new GOME-2 satellite product of 
Theys et al. (2011). Considerable improvement is achieved 
by including in GEOS-Chem the HOBr+HBr reaction in 
aerosols. This brings model BrO into the observed range 
and reproduces major features of the satellite observations 
including the increase with latitude and the seasonal max- 
ima in northern mid-latitudes in winter and in the Arctic in 
spring. Laboratory studies provide ample evidence for the 
HOBr+HBr reaction but offer limited quantitative informa- 
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tion to constrain models. The reactive uptake probability ap- 
proach used here (y = 0.2 for sulfate and sea-salt aerosols, 
y =0.1 for ice crystals) is probably too simplistic. There 
is a need to better characterize this reaction over a range of 
aerosol types and temperatures relevant to the troposphere. 

Tropospheric bromine chemistry in GEOS-Chem de- 
creases the present-day global burden of tropospheric ozone 
by 6.5 %. This is due in part to loss of NO x by BrN03 for- 
mation and hydrolysis, and in part to catalytic loss of ozone 
by HOBr formation and photolysis. Ozone mixing ratios de- 
crease by < 1-8 nmol mol -1 depending on region and sea- 
son, with the largest effects in the extratropical Northern 
Hemisphere in spring. The global mean OH concentration 
decreases by 4 %. These changes in ozone and OH appear to 
ameliorate previous GEOS-Chem biases in comparisons to 
observations, although more work is needed to evaluate the 
simulation of ozone. 

We find that bromine chemistry significantly improves 
the ability of global models to reproduce surface ozone ob- 
servations from the turn of the 20th century. BrO concen- 
trations in GEOS-Chem are comparable in the present-day 
and pre-industrial atmospheres, as the effect of decreased 
methane (suppressing the BrO x sink from the Q-LO + Br 
reaction) is compensated by the effect of decreased aerosol 
(suppressing the BrO x source from the HBr+HOBr re- 
action). Comparison of pre -industrial simulations with the 
GEOS-Chem and p-TOMCAT CTMs to 1876-1886 obser- 
vations at Montsouris (France) show mean overestimates 
of 6 nmol mor 1 in GEOS-Chem and 3 nmol mol -1 in p- 
TOMCAT without bromine chemistry. Including bromine 
chemistry has similar effects in both models, fully correcting 
the bias relative to observations in p-TOMCAT and reduc- 
ing it by half in GEOS-Chem. It suppresses the winter-spring 
model maximum, better representing the aseasonal behavior 
of the observations. Bromine chemistry in GEOS-Chem de- 
creases the natural surface ozone background over the US by 
3 nmol mol -1 on an annual mean basis, a significant incre- 
ment when assessing the effect of anthropogenic emissions 
on ozone exposure. Because the effect of bromine chemistry 
on ozone is similar in the present-day and pre-industrial at- 
mospheres, we find that the global radiative forcing from an- 
thropogenic tropospheric ozone is negligibly affected. How- 
ever, the impact on radiative forcing may be sensitive to 
changes in sea salt debromination and assumptions about het- 
erogeneous chemistry that require further investigation. 

In contrast to the similarity between BrO concentrations 
in the pre-industrial troposphere and present-day, we find 
that Br atom concentrations are 40 % higher (global mean) 
in the pre -industrial troposphere because of lower ozone. 
It has been hypothesized that Br could provide the main 
atmospheric oxidant for conversion of elemental mercury 
Hg(0) to Hg(II), with a linear-to-quadratic dependence on 
Br concentrations (Holmes et al., 2006, 2010). If so, anthro- 
pogenic ozone would have an important impact on the at- 
mospheric lifetime of mercury. Based on best estimates for 
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Hg-Br kinetics, we find in GEOS-Chem that the lifetime of 
Hg(0) against oxidation to Hg(II) by Br atoms was 4 months 
in the pre-industrial atmosphere as compared to 7 months at 
present. Most of the anthropogenic increase in tropospheric 
ozone has taken place since 1950, and background surface 
ozone concentrations at northern mid-latitudes have contin- 
ued to increase in recent decades. Thus historical mercury 
emissions from human activity (mining, industry, combus- 
tion) may not have dispersed globally to the extent previ- 
ously thought, but instead deposited mostly to the northern 
mid-latitudes. This could possibly explain observations of 
elevated mercury in North Atlantic subsurface waters (So- 
erensen et al., 2010) and the recent trend of decreasing at- 
mospheric concentrations over the North Atlantic (Slemr et 
al., 2011). It should be emphasized, however, that the role 
of Br as a global mercury oxidant remains hypothetical. Our 
results stress the need for further research to better establish 
this role. 
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